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Abstract — This paper studies a novel underlay MIMO cognitive 
radio (CR) system, where the instantaneous or statistical channel 
state information (CSI) of the interfering channels to the primary 
receivers (PRs) is completely unknown to the CR. For the single 
underlay receiver scenario, we assume a minimum information 
rate must be guaranteed on the CR main channel whose CSI is 
known at the CR transmitter. We first show that low-rank CR 
interference is preferable for improving the throughput of the 
PRs compared with spreading less power over more transmit 
dimensions. Based on this observation, we then propose a rank 
minimization CR transmission strategy assuming a minimum 
information rate must be guaranteed on the CR main channel. 
We propose a simple solution referred to as frugal waterflUing 
(FWF) that uses the least amount of power required to achieve 
the rate constraint with a minimum-rank transmit covariance 
matrix. We also present two heuristic approaches that have been 
used in prior work to transform rank minimization problems 
into convex optimization problems. The proposed schemes are 
then generalized to an underlay MIMO CR downlink network 
with multiple receivers. Finally, a theoretical analysis of the 
interference temperature and leakage rate outage probabilities at 
the PR is presented for Rayleigh fading channels. We demonstrate 
that the direct FWF solution leads to higher PR throughput even 
though it has higher interference "temperature (IT) compared 
with the heuristic methods and classic waterflUing, which calls 
into question the use of IT as a metric for CR interference. 

Index Terms — MIMO cognitive radio, spectrum underlay, 
outage probability, interference mitigation. 



I. INTRODUCTION 

Cognitive radios (CRs) have gained prominence as an 
efficient method of improving spectrum utilization by allowing 
coexistence with licensed networks, which is denoted as dy- 
namic spectrum access (DSA). One popular variant of DSA is 
known as spectrum underlay, where underlay cognitive trans- 
mitters (UCTs) operate simultaneously with licensed (primary) 
users, but adapt their transmission parameters so as to confine 
the interference perceived at the primary receivers (PRs) to a 
pre-specified threshold 1 1 1. Therefore, a fundamental challenge 
for the CR is to balance between maximizing its own transmit 
rate and minimizing the interference it causes to the PRs. 
In CR networks with single-antenna nodes, this is usually 
achieved by exploiting some knowledge of the interfering 
cross-channels to the PRs at the UCT and performing some 
admission control algorithms with power control |j2]. 

If UCTs are equipped with multiple antennas, the available 
spatial degrees of freedom can be used to mitigate interference 
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to the PRs during transmission to the underlay receivers. 
Multi-antenna CR networks have recently received extensive 
attention, assuming some knowledge of the interfering cross- 
channels to the PRs at the UCT, either perfect PR cross- 
channel state information (CSI) L3J-[5J, perturbed PR CSI [6|- 
[8 |, or statistical PR CSI S-lUl. However, the UCT may not 
have the luxury of knowing the CSI of the cross links to the 
PRs, as the primary system would not deliberately coordinate 
the collection of CSI for the CR system. 

In this work, we consider the novel scenario where both 
the realizations and distribution of the PR cross-channels 
are completely unknown at the CR, thereby precluding the 
overwhelming majority of existing spectrum underlay schemes 
in the literature [|3], El- Such a scenario of completely- 
unknown PR CSI is relevant in a number of instances, for 
example, when the PR transmits intermittently and therefore 
stymies attempts to learn the cross-channel, when channels are 
varying rapidly over time, when the PT and PR do not employ 
time-division duplexing as assumed in fW\, fill among others, 
or when there are a plurality of active PTs/UCTs and it is 
impossible to indirectly estimate specific channels. 

Specifically, we propose a rank minimization transmission 
strategy for the UCT while maintaining a minimum informa- 
tion rate on the CR link, and we present a simple solution 
referred to as frugal waterflUing (FWF) that uses the least 
amount of power required to achieve the rate constraint with 
a minimum-rank covariance matrix. In the context of MIMO 
interference channels (for which the CR underlay network 
is a special case), rank-minimization has been shown to 
be a reinterpretation of interference alignment llT4l . but this 
approach requires knowledge of interfering cross-channels and 
treats the overall system sum rate or degrees-of-freedom as 
the performance metric, assumptions which are both markedly 
different from the underlay CR scenario we consider 

We also describe two heuristic approaches that have been 
used in prior work to transform rank minimization problems 
(RMP) into problems that can be solved via convex optimiza- 
tion. These approaches approximate the rank objective func- 
tion with two relaxations, one based on the nuclear norm [22], 
and the other on a log-determinant function |23|. We show 
theoretically and via numerical simulation that minimizing the 
rank of the UCT spatial covariance matrix leads to the highest 
PR throughput in general Rayleigh-fading channels, compared 
with spreading the transmit power over more dimensions. 
Furthermore, our simulations indicate that FWF provides a 
higher PR throughput than the nuclear-norm and log-det 
heuristic solutions, even though FWF has a higher interference 
"temperature" (IT). This suggests that the commonly used 
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IT metric does not accurately capture the impact of the CR 
interference on PR performance. Instead, we propose a metric 
based on interference leakage (IL) rate that more accurately 
reflects the influence of the CR interference. 

This paper is organized as follows. The underlay system 
model is introduced in Section PR CSI-unaware UCT trans- 
mit strategies for a single UCR are presented in Section |lll] 
The generalization to the underlay downlink with multiple 
UCRs is shown in Section [V] A random matrix-theoretic 
analysis of the primary outage probability due to the proposed 
strategies is given in Section |VT] The penultimate Section IVIII 
presents numerical simulations for various underlay scenarios, 
and we conclude in Section [Villi 

Notation: We will use CA/'(0, Z) to denote a circularly 
symmetric complex Gaussian distribution with zero mean and 
covariance matrix Z, £{•} to denote expectation, vec(-) the 
matrix column stacking operator, (•)^ the transpose, (•)^ the 
Hermitian transpose, (•)^^ the matrix inverse, Tr(-) the trace 
operator, |-| or det the matrix determinant, diag(a) a diagonal 
matrix with the elements of a on the main diagonal, |A| j j the 
element of A, T{x) the gamma function, and I is the 
identity matrix. 

II. System Model 
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Fig. 1 . Cognitive radio network with a multi-antenna underlay CR transmitter, 
K underlay receivers, and a single MIMO PT-PR pair 



A generic MIMO underlay CR network with K multi- 
antenna underlay receivers is shown in Fig. 1, where a primary 
system and an underlay CR system share the same spectral 
band. Since the UCT transmit strategies are independent of 
the cross-channels to primary users, the numbers of PRs and 
PTs and their array sizes can be made arbitrary; however to 
simplify notation we will consider a solitary multi-antenna 
PT-PR pair To introduce the problem we first consider the 
scenario with a single UCR, and generalize to the case of 
K>lm Sec. [V] 

We consider a multi-antenna UCT equipped with Na an- 



A^s-antenna underlay cognitive receiver (UCR). The PT is 
equipped with Np antennas and transmits signal Xp £ C^p ^ ^ 
to the A^r-antenna PR. Thus, the UCR observes 



GiXs + G2Xp + xig 



(1) 



where Gi G C'^=^'^% G2 G C'^^^'^" are the complex MIMO 
channels from the UCT and PT, and ^ CA/^(0,a-^I) is 
complex additive white Gaussian noise. We assume Gaus- 
sian signaling with zero mean and second-order statistics 
£{xsX^} = Qs, and the average UCT transmit power is 
assumed to be bounded: 

Tr(Q,) < Ps. 



The signal at the PR is given by 



(2) 



Yp = HiXp + H2XS + Hp, 

where Hi e C^^^^-.Ha £ C^-^^" are the channels from 
the PU and CR transmitters (assumed to be full-rank), and 
Hp ^ CA/'(0, cTpI) is complex additive white Gaussian noise. 
The primary signal is also modeled as a zero-mean complex 
Gaussian signal with covariance matrix Qp and average power 
constraint Tr(Qp) < Pp. We will assume Qp is fixed and the 
channels are mutually independent and each composed of i.i.d. 
zero-mean circularly symmetric complex Gaussian entries, and 
focus our attention on the design of the UCT transmit signal. 

We assume there is no cooperation between the PT and UCT 
during transmission, and that both receivers treat interfering 
signals as noise. The network is essentially an asymmetric 2- 
user MIMO interference channel, where the UCT attempts to 
minimize the interference to the PR, but no such reciprocal 
gesture is made by the PT. The interference covariance matrix 
at the PR is 

Kp = H2Q.Hf. (3) 
Define the interference temperature at the PR as ll3l- lfT2l 



Tp(Q.) -Tr(Kp 



(4) 



tennas, which transmits a signal vector x^ e 
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Without knowledge of H2 or its distribution, the UCT cannot 
directly optimize the PR interference temperature or outage 
probability as in existing underlay proposals ll3l- lfT2l . To 
our best knowledge, precoding strategies and performance 
analyses for MIMO underlay systems with completely un- 
known primary CSI have not been presented in the literature 
thus far In addition to |3|-|13| not being applicable, the 
blind interference alignment method for the 2-user MIMO 
interference channel lITSl is also precluded since it requires 
knowledge of the cross-channel coherence intervals, which we 
assume is also unknown. In lfT3l . a blind underlay precoding 
scheme is proposed where the MIMO CR iteratively updates 
its spatial covariance by observing the transmit power of a 
solitary PT. The UCT attempts to infer the least-harmful spatial 
orientation towards the PR, but requires that the PT employ 
a power control scheme monotonic in the interference caused 
by the CR, and that the cross-channel remains constant during 
the learning process. In contrast, we investigate simple non- 
iterative CR precoding strategies which do not impose any 
restrictions on the PT transmission strategy or number of PTs, 
or cross-channel coherence intervals. 
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The PT achieves the following rate on its link: 
i?p(Q,) = log2|l + HiQpHf (Kp + a^i)-! 

Similarly, the achievable rate on the CR link is 

(Q,,) = log2 I + GiQ.Gf (K, + a^jy' 



(5) 



(6) 



where K., = G2QpG|'^ represents the interference from the 
PT. 

III. A Rank Minimization strategy for 
CSI-Unaware Underlay Transmission 

We now expound on the fundamental motivation underlying 
the UCT transmission strategies proposed in this work. As we 
have seen, due to a lack of knowledge of H2 or its distri- 
bution, the UCT cannot directly optimize the PU interference 
temperature. Hence, we propose an alternative transmission 
strategy where the UCT tries to minimize a measure of the 
interference caused to the PR in a "best-effort" sense, while 
achieving a target data rate to the UCR. Assuming that Qp is 
fixed, we first show that in the clairvoyant case where the UCT 
has some knowledge of the channel to the PR (H2), a rank- 
1 UCT covariance matrix Qs causes least interference to the 
primary link, which is described in the following proposition. 

Proposition 1: The optimal solution to the clairvoyant 
problem 



max £ {Rp{Qs)} 



S.t. 



Rs (Qs) = Rb 
Tr(Q,) < Ps 
Q. h 0. 



(7a) 

(7b) 
(7c) 
(7d) 



that maximizes the average PR rate is of rank one, i.e., 
rank(Q*) = 1. 

Proof: Since Qs ^ 0, it can be expressed as = 
UAU^, where A is the diagonal matrix of eigenvalues of 
Qs and U is the unitary matrix with columns consisting of 
the eigenvectors of Qs. Defining H2 = H2U, it follow from 
Lemma 5 in [16] that the distribution of H2 is the same as 
that of H2. As a result, the average PU rate can be expressed 
as 

i?p(Qs) = $(A) 



l0g2 



dct ( I + HiQpHf (62 AHf 



all 



Thus, the problem we considered is essentially equivalent 
to constructing the diagonal matrix A with real nonnegative 
entries so as to maximize $(A) under the constraint Tr(A) = 

From 1251 Lemma 3],|l26l, we have that $(A) is a convex 
function of A. Note that given any permutation matrix 11, we 
see (using Lemma 5 in |fT6l again) that 

$(nAn^) = $(A). 

From convexity, we have 



\ n / n 



where we have used Jensen's inequality. From the transmit 
power constraint, we have ^j-j nAII-^ = {Ps/Na)lNa- 
Thus, we have proved that the least PU rate is obtained by 
A = {Ps/Na)lNa- Further, due to convexity, we can argue 
that the largest PU rate is obtain by a point farthest away from 
A = {Ps/Na)lN^. Thus, we want A* = diag(A5;, . . . , A^J 
that satisfy ftli 



max 



Na 



s. t. ^ = Ps 



Now, 




where we used (p^) — 7^" = 1^ ^nd the equality 

is satisfied by any (A|,...,A^ ) with all zeros except for 
one nonzero entry. Hence, we conclude that rank(Q*) — 
rank(A*) = 1. ■ 
Therefore, in the clairvoyant case where the UCT has some 
knowledge of the primary CSI, a rank-1 Qs causes least 
interference to the primary link and full-rank Qs causes most 
interferenc^H Of course, the optimal Qs will depend on the PR 
CSI, which we have assumed is unavailable. Still, the result 
motivates the use of a low-rank transmit covariance at the 
UCT. It is evident that the UCT does not actually require 
knowledge of the PR CSI to minimize the rank of Qs needed 
to achieve a rate target Rf, on the CR link. To exploit this 
observation, we henceforth pose the UCT precoder design 
problem when the PR CSI is completely unknown as 



(PO) : min rank(Qs) 

s.t. Rs (Qs) = Rb 
TV(Qs) < Ps 
Qs h 0. 



(8a) 
(8b) 
(8c) 
(8d) 



This is a rank-minimization problem (RMP), and in general is 
computationally hard to solve since the rank function is quasi- 
concave and not convex. As explained below, however, in this 
case a simple waterfilling solution can be obtained. In [24J, 
the mutual information of the CR link is maximized subject 
to an interference temperature constraint and arbitrary transmit 
covariance rank constraints, which implies knowledge of PR 
CSI and thus differs from this work. 

'This notion lias been echoed in prior art on MIMO interference channels 

El, m 
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IV. Solutions for the Rank Minimization Design 

Notice that the design problem (PO) is ill-posed in the 
sense that there are potentially an infinite number of solutions. 
Suppose that we find one minimum-rank solution to (PO) such 
that Qs = UAU^ for some unitary matrix U and diagonal 
matrix A satisfies Rs{Qs) = Rb, Tr(A) < Ps- If the required 
power Tr(A) is strictly smaller than P, , then we could find 
an infinite number of solutions by making small perturbations 
to U, which while leading to a higher power requirement, 
still would require less power than P,. Obviously, the solution 
with least power is more desirable for our underlay CR system 
in order to minimize the interference caused to the PR, and 
this solution can easily be found using the Frugal Waterfilling 
(FWF) approach described next. 

A. FWF Approach 

The FWF solution seeks to find the least amount of power 
required to achieve the CR rate target of Rb with the minimum 
rank transmit co variance Q^. The optimization problem can be 
solved using a combination of the classic waterfilling (CWF) 
algorithm and a simple bisection line search. The description 
for FWF is outlined as Algorithm IV-A. 1 below. In brief, FWF 
cycles through the possible number of transmit dimensions 
in ascending order starting with a rank-one Q^, and at each 
step computes the transmit power required to meet the rate 
constraint Ri, based on CWF. This requires a simple line search 
over the transmit power for each step. Once a solution is 
found that satisfies the transmit power constraint, the algorithm 
terminates. If no feasible solution is found for all Na transmit 
dimensions, the CR link will be in outage. 

Algorithm IV-A.l Frugal Waterfilling for UCT Rank/Power 

Tradeoff (El 

Require: Ps > 0,Rb > 

set r = rank(Gi) 

for A/ = 1 to r do 
Solve: 



t- log2 



p{M) = minTr(Qs) 
I + GiQ,Gf (K, + a2l 



Rb 



end for 

if p(r) > Ps then 
Declare outage 
else 

CWF solution : N = argmmp{M); 
FWF solution : = arq mm M; 

Qs determined by waterfilling p{N) over N largest 
singular values of (K^ + cr'^i) ^^^Gi 
end if 



The FWF algorithm was presented in brief without analysis 
by the authors in [19|, and through simulation were shown to 
be an effective transmission strategy in conventional downlink, 
wiretap, and underlay networks. While FWF finds an efficient 
solution to (PO), in general rank minimization problems are 



(9) 



difficult to solve and often require exponential-time complex- 
ity. Consequently, heuristic approximations to the matrix rank 
have been proposed as alternatives in order to yield simpler 
optimization problems. In particular, the nuclear norm [l22l and 
log-determinant |23| heuristics have been proposed in order to 
convexify RMP problems like (PO) and provide approximate 
solutions with polynomial-time complexity. In the discussion 
below, we show how these approximations can be applied to 
the RMP we consider in this paper 

B. Nuclear Norm and Log-det Heuristic 

The nuclear norm heuristic is based on the fact the nuclear 
norm (sum of the singular values of a matrix) is the convex 
envelope of the rank function on the unit ball. When the matrix 
is positive semidefinite, the nuclear norm is the same as the 
trace function. As a result, the design problem (PO) can be 
formulated as follows: 

(PI): min Tr(Qs) 

s.t. Rs (Qs) — Rb 

Tr(Qs) < Ps 

Qs h . 

The nuclear norm heuristic (PI) is a convex optimization 
problem and can be solved using the CWF algorithm together 
with a bisection line search (similar to FWF). It is well known 
that under the CWF algorithm, the lowest transmit power is 
achieved when rank(Qs) is chosen as large as possible (up to 
rank(Gi)). This is clearly contrary to the rank-minimization 
design formulation, which indicates that the nuclear norm 
approach for this problem is a poor approximation. 

Using the function logdet(Qs + 51) as a smooth surrogate 
for rank(Qs), the log-det heurisitic can be described as 
follows: 

(P2) : min log det (Qs + 51) 
s.t. Rs (Qs) = Rb 
Tr (Qs) < Ps 

Qs h 0, 

where 5 > can be interpreted as a small regularization 
constant (we choose 5 = 10~^ for numerical examples). 
Since the surrogate function logdet(Qs + 51) is smooth on 
the positive definite cone, it can be minimized using a local 
minimization method. We use iterative linearization to find 
a local minimum to the optimization problem (P2) ll23l . Let 
Qs denote the fcth iteration of the optimization variable Qs. 
The first-order Taylor series expansion of logdet(Qs + 51) 
about Qi'^'' is given by 

logdet(Qs + <5I) « logdet(Q('=) +(51)+ 

Tr[(QW+JI)-i(Qs-QW)]. 

Hence, we could attempt to minimize logdct(Qs + 51) by 
iteratively minimizing the local linearization (fTTT i. This leads 
to 

Qi'^+i) = argmin Tr[(QW + (5I)-iQs] . (12) 

If we choose qI"^ ~ I, the first iteration of (fT2t is equivalent 
to minimizing the trace of Qs. Therefore, this heuristic can 



(10) 



(11) 
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be viewed as a refinement of the nuclear norm heuristic. As a 
result, we always pick Qf^ = I, so that Qi^-* is the result of 
the trace heuristic, and the iterations that follow try to reduce 
the rank of qI^-* further 

Note that at each iteration we will solve a weighted trace 
minimization problem, which is equivalent to the following 
optimization problem 

(P2-1) : min Tr (F^AF) 

s.t. log2 det (I + F^RF) = i?6 (13) 
Tr(F^F) < Ps. 

where A = (qI'^^ R = Gf (GaQpGl^ +a2l)-iGi. 
This is a Schur-concave optimization problem with multi- 
ple trace/log-det constraints. From Theorem 1 in [24|, the 
optimal solution to problem ( fT3T l is F* = A^^/^US, and 
Q(fe+i) ^ -p*F*" is an optimal solution to (O, where 
A-i/2 = UaA^^^'^U^, Ua and Aa are defined in the 
eigen-decomposition A = UaAaU^, U is a unitary matrix, 
and S = diag(y^) is a rectangular diagonal matrix. 

Substituting the optimal solution structure F* into (fTsT l. we 
have the following equivalent problem 

(P2-2) : min Tr (SS^) 

s.t. log2 det (l + S^U^RUS^ = i?6 (14) 
Tr (U^A^^USS^) < Ps 

where R = A~^/-^RA~^/^. It is found that the equivalent 
problem ( fT4] i is essentially equivalent to the converse formu- 
lation 

(P2-3) : max log2 det (l + S^U^RUS^ 

s.t. Tr(SS^)=Po (15) 
Tr (U^A-^USS^) < P^. 

This is because both formulations ( fT4] i and (fTsT i describe the 
same tradeoff curve of performance versus power Therefore, 
the quality-constrained problem (P2-2) can be numerically 
solved by iteratively solving the power-constrained problem 
(P2-3), combined with the bisection method. 

The problem formulation in (P2-3) is a Schur-convex opti- 
mization problem with two trace constraints. Using Theorem 
1 in 1.24] again, if we let R — Up^A^^U^ denote the eigen- 
decomposition of R, then the optimal unitary matrix U will 
be chosen as Up^. Denoting a ~ diag(U^A^^U) and letting 
1 > 2 — ' ' ' — ^R jv represent the diagonal elements 
of Af^, the optimal power allocation can be shown to have the 
form of a multilevel waterfilling solution: 

^] , ^ = l,...,iVa (16) 

where is the ith element of a, and /i, v can be shown to be 
the nonnegative Lagrange multipliers associated with the two 
power constraints. The algorithmic description for the log -det 
heuristic approach is outlined in Algorithm IIV-B.II 



Algorithm IV-B.l Iterative log-det heuristic Algorithm for 
rank-minimization problem 
Require: Ps > 0,Rb> 0, 

set 6 = 10-6, A = 10-3, fc = 0, 

Qi") =1, R=Gf(K,+a2l)-iGi. 

repeat 

A = (Qi'^' + 5I)-\ R = A-1/2RA-1/2, 

eig(R) = Uj^Aj^Ug. 

set U = Uj^, a = diag(U^A-iU). 

Solve: 

min l^p 

S. t. log2 (^{1 + P^XR.,?j ^Rb 

a^p < Ps. 

S = diag(p) 
F = A-1/2US 

Q(fe+1) ^ ppff 

until log2(det(Qi''' +(5I))-log2(det(Qi''+^' < A 

if a^p > Ps then 
Declare outage 
else 

p = a^p/P,; Q, = Qi'+'^ 
end if 



V. Underlay CR Downlink 

In this section we extend the blind underlay precoding 
paradigm to a MIMO underlay downlink network with K 
UCRs. We consider a modified block-diagonalization precod- 
ing strategy |27| where multiple data streams are transmitted 
to each UCR. Let each UCR be equipped with Ns antennas for 
simplicity, although the proposed precoding schemes hold for 
heterogeneous receiver array sizes as long as the total number 
of receive antennas does not exceed Na- The extension to 
the case where the UCT serves Na spatial streams regardless 
of the total number of receive antennas can be made using 
the coordinated beamforming approach ||271 . for example. The 
received signal at UCR k is now 

Yk = Gfca WfcS„,fe + ^ Gfc,iWjSuj + Gfc,2Sp + rife (17) 

where Gfe,i e C^=^^" is the main channel, Wfe e C^^xi^ 
is the precoding matrix applied to signal s^ fe G 0k for 
user k, Sp is the PT signal received over interfering channel 
Gfc,2 e C^-x^p, and ~ C7V(0,crp) is additive Gaussian 
noise. The UCT transmit covariance per UCR is now Qk,s = 
WfcW^, and the overall UCT transmit covariance assuming 
independent messages is Qg = X]fc=i Qfe s- 

We assume each UCR has a desired information rate of i?fe, 
and adopt the "BD for power control" approach in [27 Sec. II- 

1 /2 

B]. Letting Wfe = TfcA^ , it is possible to separately design 
the beamforming matrix Tfe and diagonal power allocation 
matrix Afe per user to achieve rate Rk in a two-step process. 
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Let 



G_, = [ Gi,i 



G 



fc-1,1 



'fc+1,1 



G 



K,l 



represent the the overall UCR downlink channel excluding 
the fc*'' user First, a closed-form solution for the unit-power 
beamforming matrix of user k is obtained from the 
nullspace of G_fc. To achieve this, from the SVD G_fe = 
U_fcS-fe[ V-k,! V-k.o the last {Na - h) right sin- 
gular vectors contained in V„fc^o can be used to construct Tfc 
||271 . The BD strategy therefore completely eliminates intra- 
UCR interference on the underlay downlink, and the residual 
interference-plus-noise covariance matrix at UCR k is 



Zfc — Gfc^2QpG^2 + ^sl- 



(18) 



Proceeding to the power allocation step, let 
rank(Gfc_iTfc) — Vk for user fc's effective channel, and 
assume Ik — rk- Consider the SVD of user fc's pre-whitened 
effective channel 

Z-^/'Gfe,iTfe = UfcAfcVf 

where Afc = diag (Afc^i, . . . , Afe ) is the power allocation 
matrix. While ||27l computes using the classic waterfilling 
algorithm in order to minimize the power required to achieve 
rate Rk, we can instead apply any of the other schemes 
discussed in Sec. Hill such as FWF. Due to the subadditivity of 
the rank function, reducing the rank of the per-user transmit 
covariances via FWF effectively reduces the rank of the overall 
UCT transmit covariance Q^, which in turn mitigates the 
interference caused to the PR according to Proposition 1. 

VI. Primary Outage Probability 

In this section we characterize the impact of the classic 
and frugal waterfilling methods on the primary receiver perfor- 
mance assuming independent Rayleigh fading on all channels. 
Herein, the channels are mutually independent and are each 
composed of i.i.d. zero-mean circularly symmetric complex 
Gaussian entries, i.e., vec(Gi) ~ CAf{0,T), and the same 
distribution holds for Hi,H2, and G2. 

A first approach would be to directly analyze the PR rate 
outage probability = Pr {Rp (Qs) < T) for a target rate T, 
which is equivalent to 



Pr 



(log2 



I + HiQpHf (K 



< T 



(19) 



where Qs and Qp are obtained via one of the waterfilling 
methods on their noise-prewhitened channels and are there- 
fore functions of random matrices {Gi,G2} and {Hi,H2}, 
respectively. Unfortunately, the computation of ([19) is pro- 
hibitively complex since it is a non-linear function of the 
eigenvalues of four complex Gaussian random matrices (even 
if the PT applies uniform power allocation instead), and is an 
open problem to our best knowledge. Previous studies on the 
statistical distribution of MIMO capacity under interference 
usually circumvent this difficulty by assuming the MIMO 
transmitter and interferer adopt uniform or deterministic power 
allocation ll28l . ll29l . 137 j . which reduces the problem to 
one involving two complex Gaussian random matrices. As 
such, we are not aware of prior work on the statistics of 



MIMO capacity under interference where either one or both 
transmitters employ waterfilling as in our model. 

In light of the above, it is of interest to develop more 
tractable PR performance measures. One such candidate is the 
PR interference temperature outage probability (ITOP), which 
is the probability that Tp (Qs) [cf. dUi] exceeds a threshold ry: 



/p(Q„ry) =Pr(Tr(H2Q,Hf) > 



(20) 



The ITOP is appealing since the interference temperature met- 
ric is widely used in underlay systems, and can be considered 
to be the MIMO counterpart of efforts to characterize the 
statistical distribution of aggregate UCT interference in single- 
antenna networks as in ll30l . Paradoxically, however, it is seen 
in Sec. IVIII that FWF causes the highest ITOP, even though 
the average primary rate is the highest and the PR rate outage 
is the lowest when is computed using FWF. Therefore, 
a more accurate surrogate for the PR rate outage Ir is the 
interference leakage-rate outage probability (ILOP), defined 
as 

/i(Qs,77) = Pr(log2|a2l + H2Q,Hf| > v) , (21) 



and it is verified in Sec. IVIII that UCT transmission schemes 
with the lowest ILOP also minimize 7^. This is because the 
leakage rate has a direct impact on the PR rate: Rp{Qs) in 
(|5]l can be rewritten as 



i?p(Q.) = log2|a2l + HiQpHf +Kp 



-iog2 kpi- 



(22) 



where the first term is the sum rate of the virtual PT/UCT mul- 
tiple access channel (MAC) with optimal successive detection, 
and the second term is the leakage rate from the UCT. For the 
worst-case scenario where the PT is decoded first in the virtual 
MAC, decreasing the leakage rate improves the detection of 
the PT signal in the first term and simultaneously reduces the 
second term, thereby decreasing Ir- On the other hand, the 
link between interference temperature and PR rate is more 
tenuous. 

Assume the UCT transmit covariance matrix is of rank 
fc, 1 < fc < min (Na, Ng), where k is determined by the choice 
of waterfilling scheme to achieve rate Ri, over the pre-whitened 
UCT channel Gi ^ (K^ + ci^l) "^/^d. Assume GfGi 

is of rank d', with non-zero ordered eigenvalues {ailf^^. 
Waterfilling yields a diagonal Qs with entries liT6l 



n 

a; 



l,-.-,^a, 



(23) 



where the waterfilling level ^ is a function of Ps, crl, 
and a. = (ai,...,ad') ll3Tl . lf33l . Similar arguments 
hold for the underlay downlink covariance Qs designed 
for sum-rate target Rk and aggregate channel G = 



Z-^/'GiaTi 



Z^^^^Ga'^Tx 



A. Interference Temperature Outage Probability 

Noting that the ITOP and ILOP ©-(EB are still functions 
of four complex Gaussian matrices, we first develop bounds 
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on the ITOP as follows. We assume H2 is full rank such that (^5^ — 
rank(H2) ~ d — mm{Nr, Na), which holds with probability 
1 under i.i.d. Rayleigh fading. Define r = min (fc, d) and let 
Xi (A) denote the i*^ ordered eigenvalue of A in descending 
order Starting with the commutativity of the trace operator, 

Pr (Tr (Hf H2Q,) > 77) (24) 
Pr (Tr (Hf H2Q,|Q^=(p^/A,^)i) > 77) (25) 



fr'+'^F (j + k), and the normalization 



constant Ki [38, eq. (25)] and the entries [Z] 



eq. (27)] 



< 



< 



< 




PrjAi (q.)^A, (HfH2) >j? 



(26) 



(27) 



(28) 



(29) 



i=l 



where in ( |25] l we eliminate dependence on Hi by assuming 
the PT adopts uniform power allocation (Qp = {Pp/Np)T) 
such that Ks = {Pp/Np)G2G2^ which is a worst-case 
interference scenario at the UCR according to Proposition 1 
and potentially increases the power expended by the UCT; in 
(|27| i Qs is the statistical waterfilling solution where /z in (|23] i 
is a function of the statistics of {ai}f^^ and offers nearly the 
same performance as instantaneous waterfilling llSTI - lBS] ; and 
the inequality in (l28T l follows from the bound on the trace of 
a product of Hermitian matrices (34]. 

We now define the ordered eigenvalue vectors 
h = (Ai(HfH2),...,A, (H|^H2)) and q = 
^Ai ^Qs^ , . . . , Ar I^Qs^y Observe that the overall joint 
density of these random eigenvalues is given by the product 
of the individual joint densities: /h,q (h, q) — /h (h) /q (q) 
due to the independence of the associated channel matrices. 
Thus, the bound in (|29]l can be rewritten as 



Ip < E\ 



4 



l-F, 



<ii 



V, 



.)} 



(30) 



where Fq-^ (x) is the cumulative distribution function (cdf) of 
the largest eigenvalue qi. From ( |23] l we obtain Fq-^ (x) — 
Fai (fs/(M ~ ^)) where jj is the statistical water level. The 
cdf of ai, the largest eigenvalue of G^Gi, is given next for 
the scenario Ns > Na,Ns > Np. 

Lemma 1: 1381 Given complex Gaussian matrices X e 
C^'.xJVa^Y e C^'X^p, and {Np x Np) diagonal matrix 
P = diag(p, ...,/?), the cdf of the largest eigenvalue a,nax 
of the quadratic form (YPY^ + cr^i) x when > 
Na,N, > Np is 



i^) = Ki Ai {x) 



(31) 



where Ai{x) = Y{x) 



Y{x) 



X^ ^ T 

fc=0 

(^1) ° '' Inc+n^ 



p 



1, 



1 + px ^ 
,{p), i = Na + 1, 



are functions of the array dimensions independent of x. 

The cdf Fai (x) for other antenna array dimensions is of 
a similar form and can be found in ||38l . Now, in order 
to compute the expectation over h in (|30] l, we exploit the 
Gaussian distribution of H2 based on the following lemma. 

Lemma 2: |fT6l , ||35l If X is a {N^ x No) matrix with i.i.d. 
zero-mean unit-variance complex Gaussian elements, then 
W = XX^ follows a central complex Wishart distribution 
if Nr < Na, otherwise W = X^X is Wishart-distributed if 
Nr > Na- Given c = niax{Nr, Na) and m ^ min(A^r, -^a), 
the joint density of all m ordered eigenvalues of W is 



/a(Ai,...,A™) =i^^|Vi (A)|'[]. 



^ = X]-\ and K,, is a normalization constant 



1=1 \^^) 

where Vi (A) is a Vandermonde matrix with entries 
[Vi(A)], - ^^-1 
independent of A Il35l eq. 7]. 

For the case where rank(Qs) = fc is greater than 
rank(H2) = d, i.e., r = d, the term Xli=i ™ ( l30b involves 
all d ordered eigenvalues {hi, . . . , hd} and the associated joint 
density function /h (h) is given in Lemma|2] Thus ( [30l ) yields 



/p < 1 - 



where K 



K j ■■■ j |Ai(/i)||Vi(h)|[]e(/».)rfh (33) 



KiK^, h 



and the integration region is S 



{00 >hi>h2'>..->hd>Q}- This multidimensional inte- 
gral has a closed-form solution obtained from the following 
generalized Cauchy-Binet identity. 

Lemma 3: POl Lemma 2] For x ~ {xi , ■ ■ ■ , xm}, arbitrary 
integrable functions Ci(-), Ui{-), and (p(-), N x N matrix 
* (x) and M x M matrix * (x) (M < N), where $ (x) = 



Ci(x)^ C2^ 
I - l,...,N-M;j 
scalars) for i — N - 



, with entries [Ci (x) 
= l,...,iV, [C2],^. 
M + l,...,N;j 



= {xj) for 
= C2i.j (constant 
-- 1, .,7V, and 



['4' = Ui (xj), the following integral identity over do- 

main = {b > xi > X2 > ■ ■ ■ > xn > a} holds: 

y • • • y I* (x)| • I* (x)| ixk)d^ = Ml det B (34) 



fc=l 



where 

m 



/ (f (x) Cj (x) Ui (x) dx, i = 1, . . . ,N ~ AI; Vj 
C2„-, 1 = N - M + l,...,N;yj. 

(35) 

A compact solution to ( |33] | is then obtained by setting $ = 
Ai (h) with p = Pp/Np, * (x) = Vi (h), and (x) = C (^) 
in Lemma |3] 

/p < l-d!i^|Bi|; (36) 



,N„, 



[Bi 



C^{x)xy' Y{x) 



dx, i = 1, 



[Z],^^, z = 7V-M + l,.. 



.,iV-M;Vj 
(37) 
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and requires a computationally-inexpensive one-dimensional 
numerical integration of the product of elementary functions 
in 

On the other hand, when r = k < d, only the k largest 
eigenvalues are included in the term J2l=i < l30l l. which 
necessitates invoking the corresponding joint density function 
given below in Lemma |4] 

Lemma 4: f39| The joint density function of the ordered 
subset of the s largest eigenvalues of Wishart matrix W having 
m non-zero eigenvalues in total is 



/ai....,a, (A) = iCtu ^ ^ sgn (n, m) |D (A^ 



E 

n.m.s m.m.s 



I 

where n = m = {!,...,«}, each summation is a A^^-fold 
nested sum over permutations n,rj „i of the index sets as 
defined in ll39l eq. (16)] with sign determined by sqn (n, m) G 
{±1} m eq. (17)], 



(38) 



7 {Na - Nr+ + ''j.m, A^) 



(39) 



and 7(a, b) is the incomplete gamma function 

A' ■ 



Substituting 1^(2;^) = xf" ^-+"1 ™i+2 associated 



expressions into dSOl l and invoking Lemma |3] provides 
k\K 



Ip < 1 



E E 

n,cZ,fc m,d,k 



—r 1^ sgn {n,m) IB2I 



(40) 



where 

[B2^ 



J,°°^{x,)[T,{x\^ Y(x 



2J 



dx, i — 1, . 



2J 



i = d-k + l,...,N;yj. 



VII. Numerical Results 

In this section, we present some numerical examples 
to demonstrate the performance of the proposed rank- 
minimization UCT transmit covariance designs in MIMO 
cognitive radio networks. We consider MIMO cognitive radio 
networks with one primary user and one or more under- 
lay receivers. In all simulations, the channel matrices and 
background noise samples were assumed to be composed 
of independent, zero-mean Gaussian random variables with 
unit variance. In situations where the desired rate for UCT 
cannot be achieved with the given Pg, rather than indicate an 
outage, we simply assign all power to transmit signals. The 
performance is evaluated by averaging over 1000 independent 
channel reaUzations. 



A. Single Underlay Receiver Scenario 



,d — k 



(41) 

Since 7(a, b) is a standard function in MATLAB, the numer- 
ical integration of its product with two elementary functions 
in ( |4T] ) is straightforward. 

B. Interference Leakage Outage Probability 

Turning our attention to the ILOP, starting with the defini- 
tion of /; we have 

(Q„ 77) = Pr log2 i'^l + A, (QsHf H2)) > V 

(42) 

(43) 

(44) 

(45) 




5 10 1 

UCR Desired Rate (bits/s/Hz) 



5 10 15 

UCR Desired Rate (bits/s/Hz) 



« Pr(log2n^^^A,(Q,HfH2) >77 
< Pr (n^^^ A, (q,) a, (H2Hf ) > 2 



£;h 2" 



'nl. 



l/r 



where in ( |42] | we consider the interference-limited scenario 
which is of interest, and (l43T l is due to ll34l . 

The computation of (|45T l closely parallels that of (|30] l, by 
again separating the cases r ^ d < k and r = k < d, followed 
by invoking Lemmas [T]-[3] for the former and Lemmas I1I4I 
and [3] for the latter case, respectively. Therefore, the resulting 
closed-form bounds for the ILOP are of the form in (136) and 



(Unil, with h 
unchanged. 



(2"/ni 



. l/r 



and all other terms being 



Fig. 2. Power and dimension allocation versus UCR desired rate. 

We first consider the single UCR scenario, where each node 
is equipped with 6 antennas, and P, ~ 100, Pp = 10. Fig. |2] 
illustrates the average fractional power p and the average 
number of subchannels N required to achieve the UCR desired 
rates by CWF, FWF and log-det heuristic algorithms. It is 
shown that the trace and rank of the UCT transmit covariance 
matrix are two competing objectives, and any scheme 
which requires more power occupies fewer spatial dimensions. 
Among the three methods, CWF demands the largest spatial 
footprint, while the FWF scheme offers the smallest feasible 
number of transmit dimension. We should point out that 
the log-det heuristic algorithm for matrix rank minimization 
does not always provide the smallest transmit dimension, 
compared to FWF. This is because the log-det algorithm is 
an approximate heuristic and can only give a local minimum. 

The achieved average primary user data rates for all the 
methods is depicted in Fig. [3] As expected, the lower-rank 
UCT transmit convariance will cause lesser degradation on 
average to the PU communication link, thus resulting a higher 
PR rate in accordance with Proposition 1 . Compared to CWF, 
either of the proposed modified waterfilling algorithms or the 
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5 10 15 
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Fig. 3. Achieved PU rate versus UCT desired rate. 

log-det heuristic lead to more desirable PR rates, with the 
advantage of FWF being more pronounced as Rf, increases. 




5 10 15 5 10 15 



UCR Desired Rate (bits/s/Hz) UCR Desired Rate (bits/s/Hz) 

Fig. 4. Two metrics of PR Interference versus UCR desired rate. 

To obtain greater insight. Fig. |4] compares two metrics 
of interference at PR using different algorithms, where one 
is the newly-defined UCT-PR leakage rate, the other one 
is the commonly-used interference temperature. We notice 
that an interesting phenomenon: the two metrics gives the 
opposite trend. It is worth to point out that the commonly-used 
interference temperature metric does not accurately capture the 
interference impact caused by the UCT on the primary mutual 
information, while the interference leakage rate remedies this 
defect. 

For the statistical characterization of the proposed schemes, 
we exhibit the empirical complementary cdfs and select analyt- 
ical upper bounds from Sec. |Vl]of the interference temperature 
and leakage rate metrics in Fig. |5] for 6 antennas at all users 
and Ps — 200, Rh — 8,Pp = 40. An immediate observation is 
the conflicting trends of the leakage and temperature metrics: 
FWF causes a much greater interference temperature outage 



Empirical CCDF 




5 10 15 20 25 30 35 40 



Fig. 5. Empirical ccdf of interference temperature and leakage rate under 
CWF and FWR 

and much smaller leakage rate outage compared to CWF, and 
the superiority of one versus the other is not apparent. To 
resolve this dilemma, the corresponding empirical PR rate 
ccdfs are shown in Fig. |6] and it is clear that employing 
FWF leads to a very significant reduction in PR rate outage 
probability as compared to CWF. Furthermore, the interference 
temperature outage is again seen to be misleading regarding 
the true impact on the PR rate outage probability. Thus, FWF 
outperforms CWF in terms of both average PR rate and PR 
rate outage probability. 



Empiricai CCDF of PR Rate 
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Fig. 6. Empirical ccdf of PR rate under CWF and FWF. 



B. MIMO Underlay Downlink Scenario 

Next, we evaluate the performance of the proposed algo- 
rithms with the modified BD strategy of Sec. |V]for a MIMO 
underlay downlink system, where there are X = 3 UCRs, and 
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Fig. 7. Achieved PR rate versus per-UCR desired rate in MIMO underlay Fig. 9. Two metrics of PR interference versus UCT transmit power in MIMO 



downlink, Ps = 20dB, Pp = lOdB. 



underlay downlink with identical target rates Ri = R2 = R3 



Na ^ 12, Np ^ Nr = Ns = 4. Without loss of generality it 
is assumed that the desired rate targets for all UCRs are the 
same, i.e. Ri — R2 — Rs- Fig. [T] illustrates the achieved 
PU rate versus per UCT desired rate, when Pg = 100 or 
2QdB and Pp = 10 or lOdB. The benefit of minimizing the 
transmit covariance rank is seen to hold even for the multi-user 
downlink scenario. 




13 14 15 16 17 
CR Transmit Power P [dB] 



Fig. 8. Achieved PR rate versus UCT transmit power in MIMO downlink 
system with identical target rates Ri = R2 = R3 = 5bits/s/Hz. 

It is also of interest to see how the achieved PR rates under 
the various designs vary with the UCT transmit power, when 
the desired rate for each UCT is fixed. The simulation settings 
are the same as above, except that we fix Ri — R2 — R3 = 5 
and Ps e [lOdB, 20dB]. The results are shown in Fig.|8] with 
the corresponding leakage rate and interference temperature 
metrics in Fig. |9] Once again, FWF offers the optimal average 
PR rate and PR rate outage probability in the downlink 
scenario. 



VIII. Conclusion 

This paper has proposed a rank minimization precoding 
strategy for underlay MIMO CR systems with completely 
unknown primary CSI, assuming a minimum information rate 
must be guaranteed on the CR main channel. We presented a 
simple waterfilling approach can be used to find the minimum 
rank transmit covariance that achieves the desired CR rate 
with minimum power We also presented two alternatives to 
FWF that are based on convex approximations to the minrank 
criterion, one that leads to conventional waterfilling for our 
CR problem, and another based on a log-det heuristic. The 
CWF approach turns out to be a poor approximation to 
the min-rank objective, while the log-det approach provides 
performance similar to FWF, although FWF consistently leads 
to the highest throughput for the primary link. We also 
observed that reducing the inteference temperature metric is 
surprisingly not consistent with improving the PR throughput; 
in particular, FWF has the highest interference temperature of 
the algorithms studied, but also leads to the highest PR rate. 
As an alternative, we proposed an interference leakage metric 
that is a better indicator of the impact of the CR on the primary 
Unk. 
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